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Eksperymentalne badanie frakcjonowania izotopów węgła 
między HCO: i gazowym CO: 


The, experimental study of carbon isotope fractionation 
between HCO,” апа CO; 


Экспериментальное исследование фракционирования изотопов углерода 
между НСО;- и газовым СО, 


WSTĘP 


Znajomość frakcjonowania izotopów węgla między roztworem 
i fazą gazową jest bardzo ważna w badaniach izotopowych węgla 


w wodach naturalnych, w których mogą być obecne znaczne ilości 


2, нсо5 і соз. Z tego względu w literaturze znaj- 


dujemy wiele prac poświęconych obliczaniu i wyznaczaniu stałych rów- 


węgla w postaci CO 
nowagi różnych reakcji wymiany izotopowej: 


R 
Fee gaz ёо, roztwór <-> "2CO>gaz + 13а roztwór, 
` 


2 2 


13 


K 
3 roztwór == 120 gaz + н^2с65 roztwór, 


со, па2 + н12со 


к 
sa 1 
1350 gaz + EQ roztwór za 2с0 gaz + = 


2 3 СО, roztwór. 


3 


Pierwsze teoretyczne wyliczenia stałej równowagi dla wymiany 
izotopów węgla między gazowym со, 1 соз w roztworze podane zo- 


stały już w 1947 r. przez U г е уа [ zaś dla reakcji wymia- 
ny między gazowym i rozpuszczonym w wodzie со, - м г. 1965 
przez Thode i in. [2] ‚ Do chwili obecnej brak jest obliczeń 


stałej równowagi dła wymiany izotopów węgla między gazowym 595 
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i jonem HCO3, gdyż nie można określić podstawowych częstości 
drgań jonu HCO, Ten fakt przyczynił się do konieczności wyzna- 
czenia stałej równowaci к, na drodze badań eksperymentalnych, 

Frakcjonowaniem ESEE węgla między HCO> i gazowym со, 
zajmowało się wielu badaczy [3-10]. Uzyskane przez nich wyniki 
wykazują znaczne rozbieżności, dlatego też autorka przeprowadziła 
laboratoryjne badania frakcjonowania izotopów węgla we wspomnianym 
układzie. 

Oprócz wyznaczania wielkości frakcjonowania izotopów węgla 
w stanie równowagi izotopowej badano w niniejszej pracy kinetykę 


zachodzącej reakcji wymiany izotopowej. 


CZĘSC EKSPERYMENTALNA 


Celem zbadania wyróżnienia izotopowego przygotowano roztwór 
wodny NaHCO, 2 
używano roztworu o takim stężeniu, aby pH roztworu wynosiło 


oraz gazowy СО.. Do każdeso z eksperymentów 
ок. 8,5, gdyż wtedy w roztworze istnieją praktycznie jedynie jony 


HCOZ, natomiast nie ma jonów соч cząsteczkowego со, [11 | . 


Stężenie roztworu kwaśnego кн wynosiło ok, 50 gj Gazowy 
со, używany do reakcji wymiany ргерагомапу był z tego samego 
NaHCO., z którego przygotowywano roztwór wodny, a więc począł- 
kowy skład izotopowy węgla w fazie gazowej i roztworze był taki 
зат. Dwutlenek węgla preparowano z Мансс 


centowego Н,РО, рг]. 


з Za Pomocą 100 - pro- 

Do ау EM eksperymentu zastosowano szkłany układ 
składający się z dwu połączonych kranem к, zbiorników (гус. 1), 
Do dolnego zbiornika nalewano odpowiednią ilość roztworu kwaśne- 
go węglanu, układ szklany podiączano do aparatury próżniowej і 
po odpompowaniu go wprowadzano do niego gazowy cO, uzyskany 
z Мансо.. 

W układzie zachodziła wymiana izotopowa między kwaśnym 
węglanem i dwutlenkiem węgla (kran к, był otwarty, zaś K, zam= 
knięty). Cały układ termostatowano, Ponieważ w omawianym ekspery- 
mencie stosunek zawartości węgla w fazie gazowej do jego objętości 
w fazie ciekłej był Ściśle określony, wynosił R(RĆ<1), dlatego też 


podczas wymiany skład izotopowy węgła w HCO zmieniał się Row 
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-krotnie mniej niż w gazowym CO»: Z tego względu do badań wyróż- 
nienia izotopowego wystarczyło określić skład izotopowy wyjściowego 
NaHCO} oraz skład izotopowy gazowego со, pobranego znad roz- 
tworu. 


Ryc. 1, Układ szklany, 1 
w którym zachodziła re- 
akcja wymiany izotopowej 


Dla нсо* wartość 6 130” obliczano ze wzoru (otrzymanego 
w 
wprost z bilansu izotopów): 


13 13 13 
GL". бе кв (6 E 


13 
нсо; мансо, - 6 Garo; | (1) 


; W kilku przypadkach mierzono także skład izotopowy węgla w 
HCO, i stwierdzono, że wyniki były zgodne z obliczonymi wg wzoru 
(1). Próbki gazowego CO, znad roztworu pobierano ze zbiornika 2 


2 


(ryc. 1). Dla każdego punktu pomiarowego przygotowywano świeży 


roztwór i со. Pobraną próbkę gazowego со, oczyszczano 2 мо- 


dy і mierzono skład izotopowy węgla. Pomiary $ 13 przeprowadza= 
no na zmodyfikowanym spektrometrze mas МЕ1305 [13,14 |wieikość 
trakcjonowania izotopowego będzie podawana niżej w postaci Е Pe: 


C13c “ (4. 1) x 1000 (2) 


ee È (23сй2су 


we : 


©1367 


х 1000. 


218 Janina Szaran 


gdzie a. jest współczynnikiem frakcjonowania izotopowego: 


, (аса. m 


13_f12 
бе 1:7 


Dla б bliskich jedności E 5 wyraża się poprzez wartości безе; 


HCO 


228 513, wR 8 а a 
3 а... 


Błąd standardowy pomiarów © > wynosił 0,3 Ooo. 


WYNIKI BADAŃ I DYSKUSJA 
Kinetyka reakcji wymiany izotopowej 


W badaniach kinetyki reakcji wymiany izotopowej określano 
eksperymentalnie zależność сезе od czasu trwania reakcji wymia- 
ny (dla jednego czasu wymiany - 1 eksperyment). 

Przeprowadzono 3 serie eksperymentów: 

1. Przy ustalonej temperaturze 327 К zmieniano ciśnienie gazowe 
go СО. od 26,6 do 133 kPa (ciśnienie całkowite (co, + H,O) w 
układzie zmieniało się w granicach od 41 do 147 kPa). 

2. Przy ustalonym ciśnieniu p(co,)= 26,6 kPa zmieniano tem- 
peraturę od 292 do 320 KS 

3. Przy ustałonej temperaturze ok. 292 К i ciśnieniu p (со.) 
ок. 27 kPa badano kinetykę reakcji wymiany w dwu przypadkach: 

a) bez mieszania gazowego со, z roztworem kwaśnego ме- 
glanu , 

b) z intensywnym mieszaniem obu faz w momencie początkowym 
reakcji wymiany. 

W przypadku a) do zbiornika 1 (ryc. 1) nalewano odpowiednią 
ilość roztworu, następnie szklany układ odpompowywano i do zbior= 


піка 2 wprowadzano odpowiednią ilość gazowego СО. W końcu oba 


2* 
zbiorniki zostały .połączone przez obrót kranu к.. 

W przypadku b) do odpompowanego zbiornika 1 wprowadzano 
najpierw gazowy CO, a następnie zasysano do tegoż zbiornika od- 


powiednią ilość roztworu. Tym razem w momencie zasysania roztworu 
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zachodziło intensywne mieszanie się obu faz i reakcja wymiany prze- 
biegała tylko w zbiorniku 1. 
Wyniki przeprowadzonych badań są zebrane w tabelach 1, 2 


i 3 oraz przedstawione graficznie na ryc. 2, 3 i 4. 


EC |19) £.25,6 % ,T=327K 


© -ря 41 kPa 
о-р= 61 kPa 
в-рз 81 kPa 
о -p = 07 kPa 
4 -р= 147 kPa 


2 Й 6 8 10 12 1 tih) 


В ЖФ: Е к : 
Гус. 2. Zależność E 3c od czasu trwania reakcji wymiany przy 


stałej temperaturze 327.K i różnych ciśnieniach w układzie 


Zależność © 23с od czasu może być podana w postaci 


wzoru przedstawiającego równanie kinetyczne I rzędu [15] : 


См 


gdzie {се oznacza E 6 w stanie równowagi izotopowej, газ 
т = stałą czasową 

Z eksperymentalnych wyników obliczono wartości T (patrz 
tab, 1 i 2), a następnie sporządzono wykresy przedstawiające za- 
leżność: 

1) stałej czasowej od ciśnienia przy ustalonej temperaturze 
(ryc. 5), 

2) stałej czasowej od temperatury przy stałym ciśnieniu p 
(CO„) (ryc. 6.) 

г Wyniki badań Е wo jako funkcji czasu trwania wymiany w 
temp. 327 К przy różnych ciśnieniach gazowego со. (ciśnienie - 
pary wodnej w naczyniu reakcyjnym w tej temperaturze wynosiło 
ok. 14 kla) zostały zebrane w tab. 1 oraz przedstawione graficze 


nie na ryc, 2. 
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Tab. 1. Wyniki badań Е13с jako funkcji czasu przy ustalonej tempe- 


raturze 327 К i różnych ciśnieniach w układzie 


р(со,+н,о) [кра] t [h] ŻE [е] stała czasowa U [к] 


1,5 1,16 

147 3 2,47 4,6 
7 4,56 
1 1,00 
2 2,33 

107 5 4,33 3,7 
10 5,27 
17,5 5,55 
1,5 1,88 

81 3 3,62 2,9 
4,5 4,34 
Е 1,70 

61 2 3,01 2,5 
3 3,81 
0,5 2,05 
1 3,59 
2 4,53 

41 4 5,05 1,0 
7,5 5,29 


14,3 5,66 
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Eo 28,5 Yoe , plCO,) + 26,6kPa 


2 4 - 8.0 т th 


Ryc. 3. Zależność Е 132 od czasu trwania reakcji wymiany przy 
stałym ciśnieniu со, i różnych temperaturach 


W badanym układzie w chwili początkowej skład izotopowy мев- 
la w obydwu fazach był jednakowy, w miarę upływu czasu zwiększa 


się- E Sie zmierzając do wartości Ё a че odpowiadającej stanowi rów- 


nowagi izotopowej, W określonej temperaturze сетер ma oczywiście 
Ściśle określoną wartość (w temp. 327 K тата ZA %), jednakże 
szybkość osiągania stanu równowagi izotopowej jest zależna od ciś- 
nienia. Wielkością charakteryzującą tę szybkość jest ARA czasowat, 
2. За а 


0,63 Eo WE. Z dobrym przybliżeniem można przyjąć, że czas usta- 


która jest czasem, po upływie którego układ osiąga 


lania się równowagi izotopowej wynosi 5T, Jak widać z гус. 2 i 5, 


stała czasowa т rośnie ze wzrostem ciśnienia, 


EC | {ы 
n Ex 
a 
Е. = 5,5% 
T s 292K 
р(С0,) = 27kPa 
1 2 3 буть Ке 


Ryc. 4. Zależność € зе od czasu trwania reakcji wymiany; а - brak 
mieszania faz, b = intensywne mieszanie faz w początkowym stadium 


reakcji wymiany 
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Т . 327 К 


02 04 0608 10 12 K 16 18 26 22 2A 7 
Я х 10 
Ryc. 5. Zależność stałej czasowej T od cisnienia w układzie 


Ryc. 5 pokazuje, że wartość Ec 130 maleje ze wzrostem tempe 


ratury, Zależność ta będzie omówiona w dalszej części niniejszej 
pracy. Związek między stałą czasową GEAR temperaturą przy stałym 
ciśnieniu co, jest przedstawiony graficznie na гус, 6. 


р(С0;) = 26,6 kPa 


20 280 50 о 30 


Рус. 6. Zależność stałej czasowej % od temperatury 
Jak widać, dia przedziału temperatur niniejszego eksperymentu (292- 


-327 K) uzyskano liniową zależność między © i temperaturą, Pro- 
sta ta może być opisana. równaniem: 
T [һ]- 2,26 - 0,0223 (T - 273) я (6) 
Ryc, 4 ilustruje zależność zie = f(t) dla tych samych war- 
tości ciśnienią i temperatury. Krzywa a) dotyczy wymiany izotopowej 


zachodzącej bez mieszania gazowego со, i roztworu kwuśne= 
* 


— 
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л. 


go węgłanu, zaś b) - tejże wymiany, gdy w początkowym stadium wy= 
miany zachodziło silne mieszanie.obu faz. Widać wyraźnie, że w przy- 


padku mieszania faz wymiana izołopowa zachodzi znacznie szybciej 


ч 1 ‚Ж 
Tab. 2, Wyniki badań Е 3e jako funkcji czasu trwania reakcji 


wymiany izotopowej przy ciśnieniu СО, równym 26,6 kPa 


2 
i różnych temperaturach 


T [к] t Гь] €**c[ eo] stała czasowa %[h] 
1 3,29 i 
292 2 5,91 1,8 
4 7,62 
0,5 1,47 
1 Эт 
1,75 5,87 
298 4 7,13 197 
8,5 7,88 
#6 8,45 
18 8,42 
0,5 1,72 
310 з 4,09 155 
3 6,04 
. 0,5 ; 2,05 
1 3,59 
327 2 4,53 1,0 
4 5,05 
7,5 5,29 


14,3 5,66 
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Tab, 3. Wyniki badań © а јако funkcji czasu trwania reakcji 
wymiany izotopowej w przypadku bez mieszania faz a) 


i z mieszaniem faz w początkowym stadium wymiany b) 


a) bez mieszania faz b) z mieszaniem faz 
t [№] €13c [о] t [nr] ee °/со | 
3,29 0,2 4,72 
5,91 0,5 6,72 
4 7,62 0,7 7,51 
1,2 8,23 
2,1 8,45 


A 


Tab. 4. Wyniki badań Е 130 w stanie równowagi izotopowej w prze- 
dziale temperatur od 273 do 325 K 


T [x] CZE [Кое | 


273 10,10 
283,5 9,35 
293 8,34 
298 7,95 
299,5 7,90 
303 7,45 
308 7,03 
319,4 5,72 
324 5,68 


325 5,66 
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niż bez mieszania (nieznaczna różnica geometrii układu dla fazy ga- 
zowej nie ma tu większego znaczenia), Oszacowane stałe czasowe 
wynoszą; т = 1,7 h, zaś "0,2 0,25 п, 

Czyskanych wyników dotyczących kinetyki reakcji wymiany 
w badanym układzie nie można porównać z danymi literatury, po- 
nieważ w literaturze ukazała się tylko jedna praca (w r. 1940) do- 
tycząca kinetyki reakcji wymiany izotopowej w omawianym układzie 
[16 | Autorzy tej pracy podają, że reakcja wymiany izotopów węgla 
między HCO, i gazowym со, osiąga stan równowagi w ciągu kilku 
minut To oszacowanie nie może być brane pod uwagę, gdyż autorzy 
ci przyjęli błędne założenie, że równowaga izotopowa osiągnięta 
jest wtedy, gdy w нсоз i СО, jest jednakowy skład izotopowy wę- 
gla. Z prac wielu innych autorów wynika, że w badaniach ќгаксјопо- 
wania izotopów węgla w omawianym układzie, w stanie równowagi 
izotopowej, próbki były pobierane po kilkunastu godzinach е] lub 
po kilku dniach Fa 10 | ad chwili rozpoczęcia wymiany. 


Zależność temperaturowa frakcjonowania izotopów węgła w 


stanie równowagi izotopowej 


Pomiary wielkości Е e w słanie równowagi izotopowej były 
przeprowadzone w przedziale temperatur od 273 do 325 К. Układ, 
w którym zachodziła wymiana izotopowa, był termostatowany tak dłu- 
go, aż została osiągnięta równowaga izotopowa (t > 57 ), po czym 


pobierano do badań gazowy CO,. Dla badań w innej temperaturze 


przygotowywano do EEE „c świeży roztwór i со,» 
Wyniki pomiarów podane są w tabeli 4, zaś w tabeli 5 zebra- 
no wyniki badań frakcjonowania izotopów węgla między нсо5 і со, 
przeprowadzonych przez jc autorów. Matematyczne modele 
opisujące zależność między Es de i i диз podane są m.in. 
przez Mooka i in, [10] : fa эс = т: < + В, ога2 рг2е2 Deinesa 
i in, [18]: EM C = AJT"Ż + В. W tabeli 5 autor podaje także 
inny model: Ж С = AT + Ba чүс гуйч liczbowe dla tych 
modeli zostały obliczone metodą najmniejszych kwadratów na podsta- 
wie danych eksperymentalnych z cytowanych prae (patrz uwagi w 
tab. 5). Dla modeli opracowanych na podstawie wyników badań ni- 
niejszej pracy autor podaje także sumę kwadratów odchyleń, skąd 


wynika, że uzyskane przez autora dane eksperymentalne są naj- 


“иэтТАчоро Mo9ojeipemy buns ezowuzo SS -~ {х 


wa > 4 - = Ре 
685 FAAA 55 (98*0 569'5) - „ 1.07(82 70021) = 3 =: весе 
"ү = “ат ү - = 
oz? 5415'0 = SS (®г*т z9T'6T) q-T(0€S 37908) = 3 За Злам 
2835 TSYT'O = 55 (є22-1,) (всоо*о 79680'0)-(PT'Oz2T*'OT) = 3 
(x 
z м 66€-822 
4 266. 26Е [Гот] 68'cz-,_L £2v6 = 3 [Гот] "ur 1 жоо 
$ @ 
$ T8E [21] лао *zołodsm z8'92-r_LZZZOT = 3 [6] “м 1 upwa 
(0 
с ры [ГУ] чолчәр z qo -=оюдем 86'Є-г-1 000256 =3 М Coe !EŁZ 
ze {2=-1,)0'0-2'6 = 3 [9] *ыр | iosnoq 
Е Let] Кэела м 
сг рер ZÁ z цоќиер Z очотоцао м 695-962 
"хоюачзм ©MOĄZOJ ID9OLGM (52=-1,)950'0-0*6 =3 [2] "ч т чаем 
ЮЖ” ш ЖОЧУ МО Тян Е М 
*ziepuz *Tqo i f «dua; ұәү2рә2л1а 
Poan (3) "2,9 А. 
. 2.3 
Cx) MUJ ү, ст omy 


ч {әлоЧйоўот. 
{#елошмол ajus м тоб боо» — оон arzpełrm м лоЧйоўо2т GTUPMOUO(DĄYCJJ цереа ричАм "S “ат 
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lepiej. opisywane przez model 32 = AZT + в... Ten model, w 
przedziale niskich temperatur, jest także dobry ze względów ргак- 
tycznych, gdyż рее zmienia się liniowo ze zmianą temperatury. 


Тако zależność została przedstawiona na гус. 7. 


ciis) 

© 

A *- Malinin i in. (1967) 
а = Mook i in, { 1974) 

z o - niniejsza praca 

6t 

5 

4 

3 

2 

1 

0 


280 00 30 30 30 2380 40 ыб о К] 
Ryc. 7. Zależność Co 36 od temperatury: linia ciągła w zakresie 
wyższych temperatur jest prostą poprowadzoną metodą najmniej- 
szych kwadratów na podstawie danych Malinina i in. [2 zal ‚ zaś 


druga - na podstawie wyników niniejszej pracy 


Linia ciągła w zakresie wyższych temperatur prezentuje ostatnio 
wspomniany model uzyskany w oparciu o dane pracy [2] ‚ zaś 
druga prosta - w oparciu o dane niniejszej pracy. Naniesione па 
wykres punkty doświadczalne wskazują na dobrą zgodność wyników 
uzyskanych przez autora z wynikami pracy Mooka i іп. [10]. 
Wszystkie przedstawione w tab, 5 modele przewidują inwersję 
znaku Sag c w dostatecznie wysokiej temperaturze. W tabeli 5 po- 
dano temperatury inwersji, obliczone na podstawie przedstawionych 
modeli oraz zmierzone eksperymentalnie. Z porównania temp. inwersji _ 
zmierzonej i obliczonej widać, że modele matematyczne, w których wy- 
stępuje człon Яг, dają zbyt wysokie temp. inwersji. Dla pozostałych 
modeli zgodność z eksperymentem jest o wiele lepsza: T inw. mie- 
szczą się w przedziale 381- 434 К, podczas gdy zmierzone wy- 
noszą: 397 K [10] i 423.К ke]. 
Temperatura inwersji obliczona na podstawie najlepiej dopasowanego 


do wyników niniejszej pracy modelu: 


6.236 = 10,17 - 0,0896 (Т - 273) 
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wynosi 387 K, a więc jest bliska uzyskanej przez Mooka [10]: 
397 к. | 
Autor dziękuje dr 8. Hałasowi za pomoc w przeprowadzaniu 
obliczeń współczynników w modelach opracowywanych w niniejszej 
pracy. 
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SUMMARY 


Two sets of experiments on the kinetics of carbon isotope 
exchange reaction between HCOJ and gaseous со, рауе Бееп 
carried out: (a) at temperature уо! 327 Kand varyine partial pres- 
sure of со, from 27 to 133 kPa, and (b) at constant pressure ot 


со» of 27 kPa, and temperature runing from 292 to 327 К. Н has 
been demonstrated that this is a first order reaction the rate of whici. 


decreases with the pressure increase, and proportional to the tempe= 
rature. The rate of reaction hishly depends on mixing conditions be- 
tween those two phases, At T = 292 K and pressure СО 2 27 kPa, 
for example, the isotope equilibruinn may be accomplished =" 8,5 
to 1,25 hours, 

The equlibrium fractionation factor, has been measured in 
a temperature ranse from 273 to 325 К. The per mil difference 
in Roj EE between HCO, and CO, can be expressed as the fol- 


EH 


lowing function of temperature: Ẹ& 5С = 10,17 - 0,0896 (T - 273). 


РЕЗЮМЕ 


Іроведены исследования кинетики изотопного обмена углеро- 
да между НСО. и газовым 20: /а/ при температуре 227 К под дав- 
дением соь or 27 по 133 кРа и /6/ при постоянном давлении CO, 
/27 кра/ а меняющейся температуре в пределах 292 - 327 К. 


` 


130 Janina Szaran 


э„——————-—ЄїГг——————— 


Обнаружено, что изотопный обмен является реакцией первого по- 
рядка, которой скорость падает с повышением давления и растет 
линейно с повышением температуры. Показано также, что переме- 
шивание фаз резко влияет на уменьшение времени устанавливания 
изотопного равновесия, При температуре T=292 К и давлении COn 
27 KPa оно уменьшается с 8,5 до 1,25 часа, 

Измерено фактор фракционирования изотопов углерода c а 
в состоянии изотопного резновесия в пределах температур 273 - 
325 К. Установлено, что этот фактор хорошо представляется фор- 
мулой 53 = 10,17 - 0,0896 (T - 273). 


Praca finansowana z tematu MR.I.5 (Uniwersytet Warszawski). 


Złożono w Redakcji 26 IV 1983 roku. 


